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摘要 


本 论文 集中 了 一 些 关 于 一 维 随机 游 动 的 研究 和 思考 。 


首先 ， 我 们 通过 求 分 布 列 、 数 学 期 望 以 及 方差 等 数字 特征 研究 了 有 限 步 一 维 随 机 游 动 。 
然后 由 有 限 推 广 到 无 限 , 研究 了 无 限 步 一 维 随 机 游 动 。 根 据 探究 得 到 的 二 项 分 布 的 项 与 其 之 
间 的 关系 ， 利 用 棣 英 弗 方法 , 导出 了 其 概率 密度 函数 。 


随后 又 分 别 对 一 维 随机 游 动 增加 反射 壁 和 吸收 壁 两 个 条 件 。 同样 从 有 限 步 出 发 , 通过 求 
分 布 列 来 寻找 规律 。 对 比 增加 条 件 前 后 的 情况 ， 我 们 从 杨辉 三 角 的 变形 中 寻找 到 一 般 算 法 ， 
并 类 推 到 无 限 步 。 


最 后 对 全 文 进 行 了 总 结 。 


关键 词 
随机 游 动 杨辉 三 角 反射 壁 吸收 壁 
Abstract 


This thesis focuses on random walk in one dimension. 


Firstly, we studied simple random walk by finding out its distribution, mathematical 
expectation, variance and some other digital features. Then, extended to ‘infinite Step ， 
we further studied random walk. According to the relation between the items of the 
binomial distribution and random walk, we got its probability density function using 
De Moivre methods. 


‘Reflection wall’ and ‘absorbing wall’are respectively added to the simple random 
walk. We also started from ‘inite step’and found out its distribution in search of a law. 
From Pascal’s triangle, we got a general algorithm which can be extended to ‘infinite 


step’” comparing the situation before and after adding the conditions. 


At last, we drew a conclusion. 
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论 


于 
lll 


“4 于 而 行 ， 如 醉酒 者 ” 


向 的 醉 汉 不 知道 什么 时 候 分 


形象 而 又 生动 地 描述 了 随机 游 动 的 特点 ， 就 好 像 一 个 喝 得 坚 头 


误导 他 媒 


四 处 迁徙 的 种 群 ， 或 如 一 个 不 断 碰 撞 的 水 分 子 ， 都 体现 着 随机 游 动 的 特点 一 一 即 在 随机 性 
体现 着 规律 性 , 而 又 在 规律 性 中 蕴含 着 随机 性 。 而 本 文 将 着 重 研 究 人 简单 的 和 附加 条 件 的 一 
随机 游 动 ， 期 望 在 一 条 简单 的 数 轴 上 发 现 不 简单 的 规律 。 


什么 方向 行走 。 无 论 是 在 宏观 世界 ， 还 是 在 微观 世界 ， 或 如 一 


1. 1 一 维 随机 游 动 定义 


点 P 了 在 数 轴 x 上 作 随 机 游 动 ， 即 点 P 从 原点 0 出 发 , 每 次 向 x 轴 正 方向 , 或 x 轴 负 方向 


移动 1 个 单位 。 


1.2 有 限 步 一 维 随机 游 动 


问题 


求 点 P 经 过 n 次 随机 游 动 出 现在 数 轴 x 各 位 置 的 概率 , 以 及 点 P 出 现 位 置 的 数学 期 望 和 


方差 。 


分 析 


经 过 n 次 随机 游 动 ， 点 P 出 现在 位 置 s 的 可 能 情况 


由 表 1 可 得 , 经 过 n 次 随 
的 排列 


~ 
3 


机 游 动 , 点 P 出 现在 位 置 s 的 可 能 情况 的 排列 满足 “杨辉 三 角 


NS -10 | -9 -8 = -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 10 
1 1 1 
2 1 2 1 
3 1 3 3 1 
4 1 4 6 4 1 
5 1 5 10 10 5 1 
6 1 6 15 20 15 6 1 
7 1 7 21 35 35 21 7 1 
8 1 8 28 56 70 56 28 8 1 
9 1 9 36 84 126 126 84 36 9 
10 1 10 45 120 210 252 210 120 45 10 1 
表 1 


P 每 次 向 x 加 


点 了 总 共 经 过 n 次 移动 ， 令 点 P 向 x 轴 正 方向 移动 i 次， 其 
正方 向 移动 的 概率 为 p， 则 点 P 了 向 x 各 


轴 负 方向 移动 的 概率 为 (1-p)。 
那么 点 P 出 现在 位 置 (2i-n) 的 概率 P1 满 足 二 项 分 布 为 
EA 
邻 s=2i-n， 其 中 -n 三 s 万 n，s EZ， 则 点 P 出 现在 位 置 s 的 概率 Pj 为 
nts 2 一 s 
n 2 2 
ES (1—p) 
2 
解答 


点 P 出 现在 数 轴 x 各 位 置 概率 P 的 分 布 列 


PO0in, 1EZ, 并 
负 方 向 移动 (n-i) 次 ， 点 P 每 次 向 x 


E 


当 n 为 奇数 时 
xX —n —(n-2) -1 1 n-2 n 
(09-0 =p)" EA EN EN (a 1) nD (2) 
当 n 为 偶数 时 
—n =(ri-2) 了 多 0 2 n-2 n 
人 
点 P 出 现 位 置 的 数学 期 望 


民 据 数学 期 望 定义 ， 点 P 出 现 位 置 的 数学 期 望 为 


可 卫 十 S n—s 

0 oj 2 
和 n To 

根据 组 合 恒等式 m[ ”一 | 的 ， 变 形 得 


Et- i (0) pj -i- 


2np: 


根据 二 项 式 定 理 ， 化 简 整 理 得 


> 


i=0 


| 


nO—1 i 
i=1)? 


E(x)=2np-n 


点 P 出 现 位 置 的 方差 


民 据 方差 定义 ， 点 了 出现 位 置 的 方差 为 


n 


Var (x)= > (s— (2 np 一品 )2. 


人 


n++s 入 二 


n .本 2 = 
a (1— 7p) 
2 


> tone en 


R 据 组 合 恒等式 n.[ ”一 ] =m( 站 ， 变 形 得 


n 
Wn 


var->4 网 p(B > 的 汉人 


民 据 二 项 式 定 理 ， 化 简 得 


Var (x)=4n(n 一 1)p +4np -8 np +4 np 


Var (x) =4np (1-p) 


1. 3 无 限 步 一 维 随机 游 动 


考虑 p= 二 ，n 为 无 穷 大 偶数 的 情况 


今 n=2m，mE7Z， 其 中 心 项 为 


中 “ 抄 ( 打 


R 据 斯 特 林 公 式 Jim_n! = 27 ie ， 


其 中 心 项 可 表示 为 


1 


Rs 


. 2m..-2m 2m 
B (m)= 27T:2m:2n1 8 1 
m 2 


VY 2m-m ey 27010 .Ne 


即 其 中 心 项 可 表示 为 


B (m) = 2 
2m 


令 d 为 常数 ，dE2Z， 其 任意 项 为 


er 人 (2 二) i 


其 任意 项 与 中 心 项 的 比 


B(m+d) - ml:m! 
B(m) (m+ a)!:(m— aq)! 


1 

B(m+d) - V2n:m .eV/2n-m 
1 
泡 


两 边 取 对 数 可 得 
B(m+ d) 


全 


1 
E > A = AN a 
民 据 斯 特 林 公式 limn! =V2r ?.e” ,其 任意 项 与 中 心 项 的 比 可 表示 为 


| =(2 m+ 1):ln(m) (nt d 十 二 ] mo+ d) mat 二 mo- d) = 


B(m) 


为 无 穷 大 偶数 ，d 为 常数 “5 一 0 


根据 麦克 劳 林 展开 式 ln(1 + x) =x 一 + Rn(x) (对 <1)， 


(2 a 


售 去 二 阶 小 量 得 


B(m) 
HH 任意 项 可 表示 为 
-到 2 
Bmtd)=Blme "= 1 .sm 
V Tm 
即 其 任意 项 可 表示 为 
_ 2 到 
Bnrd)=_ 2 .on 
2m 


取 两 项 展开 


将 任意 k 项 相 加 得 


即 


k 2 


k 
| 2 
d=0 


oY 2m 


点 了 经 过 无 限 次 随机 游 动 出 现在 区 间 (a, b) 的 概率 为 


. 22 


nlim,,P1(a <x < 2) -| :e ”dx 


2m 
a 


令 y=_x_ ， 将 其 标准 化 得 


图 1 为 一 维 随 机 游 动 中 点 PH 


100 时 的 图 


im f(a <y<b) | -3 “dy 


像 。 


a 


图 1 


27:n 


bs 现 位 置 的 概率 密度 函数 {x)= 一 2 .6 


n 


, 当 n 从 1 取 到 


| 
n 
i 
2 


计算 机 拟 合 的 图 像 。 


图 3 


( 直 ， 当 n 从 1 取 到 100 时 | 


图 3 为 随机 游 动 中 p=0. 51 的 情况 PO 各 +0=f 0.51? “0.49? “， 当 n 从 1 
刀 


(2 十 x 


取 到 100 时 由 计算 机 拟 合 的 图 像 。 


随 着 p 的 增 大 ， 需 要 更 大 的 n 才能 更 好 的 吻合 ， 因 此 不 一 一 枚 举 图 像 。 


1.4 小 结 
由 上 可 知 ， 所 得 一 维 随机 游 动 中 点 P 出 现 位 置 的 概率 密度 函数 f(x)= 1 .。 
2 区 
较 好 的 体现 其 概率 分 布 ， 其 中 n 越 大 吻合 的 越 好 。 
1 
0.5 
-10 -5 0 5 10 
图 4 


2. 1 反射 壁 定义 


| 


点 了 在 数 轴 x 上 过 到 反射 壁 , 即 点 P 移动 一 个 单位 到 反弹 壁 后 , 便 向 相反 方向 移动 一 个 


单位 。 


2. 2 存在 反射 壁 情况 下 的 有 限 步 一 维 随机 游 动 


问题 


位 置 有 一 道 反 射 壁 ， 求 点 P 经 过 n 次 随机 游 动 出 现在 数 轴 x 各 位 置 的 概率 。 


分 析 


经 过 次 随机 游 动 ， 点 P 出 现在 位 置 s 的 可 能 情况 ( 表 2 是 r=3 的 情况 ) 


1 1 
1 2 1 
1 3 3 1 
1 4 6 6 
1 5 10 11 5 
1 6 15 21 21 
1 7 21 35 42 21 
1 8 28 57 78 84 
1 9 36 85 126 162 84 
10 45 121 220 297 330 
表 2 


由 表 2 可 得 , 经 过 n 次 随机 游 动 ， 点 P 出 现在 位 置 s 的 可 能 情况 的 排列 满足 “对 折 ” 的 
“杨辉 三 角 ” 的 排列 。 


-9 -8 -7 -6 -5 -4 -3 -2 = 0 1 2 4 5 6 7 8 


1 7 21 35+1 35+7 
1 8 28 56+1 70+8 56+28 
1 9 36 84+1 126+9 126+36 


10 45 120+1 210+10 252+45 210+120 


加 | 
1 1 
1 2 1 
1 3 3 
1 4 6 4+1 
1 5 10 10+1 
1 6 15 20+1 15+6 


对 比 表 1 和 表 2 可 得 ， 表 2 是 由 表 1 关于 r=3“ 对 折 ” 以 后 相 加 而 得 ( 见 表 3) 。 


不 妨 令 0<r<n,，rEZ， 那 么 点 P 出 现在 位 置 s 的 概率 Ps 满足 


TT 


当 s 二 2r-n 时 
P; (X=s) =P1 (X=s) 


当 2r-n 入 sf 时 

P; (X=s) =P, (X=s)+ P(X=2r-S) 
当 S> 工 时 
P: (X=s)=0 


我 们 可 以 把 反射 壁 看 作 是 一 面 镜子 , 那么 点 P 出 现在 位 置 s 的 概率 就 等 于 点 P 出 现在 位 置 s 
的 概率 与 位 置 s 关 于 反射 壁 r 对 称 位 置 s” 的 概率 之 和 。 


点 P 出 现在 数 轴 x 各 位 置 概率 P: 的 分 布 列 


当 n, r 同 为 奇偶 数 时 


—(n-2) … 2r-—n-2 2r-n 2r-n+2 


EA a Et (ran "(9):m ( n ja-m ( n 


r-2 je 
十 ntr nr 
ntr nr nt+r nr 二 
人 | 1p) 3 [rd Pp 2 We 2 | 加 外 3 
上 | \ 2 +1 2 
当 n, T 不 同 为 奇偶 数 时 
过 二 3 下 一 和 一世 
(n-2) 2 2 网 2r-n+2 


(rp oo 全 ap ( n ee n 
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2. 3 存在 反射 壁 情况 下 的 无 限 步 一 维 随 机 游 动 


民 据 1. 4 的 结论 和 2. 2 的 分 析 得 


当 x 私 7 时 
fo (x) =f1 (x) +f (2r-x) 
当 x>T 时 
f, (x) =0 
2.4 小 结 


由 上 可 知 ， 所 得 存在 一 道 反射 壁 的 情况 下 的 一 维 随 机 游 动 中 点 P 出 现 位 置 的 概率 密度 函 


好 _ (x-27)? 
Fd 2 
区 


1 
数 f (x) = 本 党 有 
0 


吻合 的 越 好 。 


X=>r 


3.1 吸收 壁 定义 


x 人 7 能 够 较 好 的 体现 其 概率 分 布 , 其 中 n 越 大 


点 p 在 数 轴 x 上 过 到 吸收 辟 ， 即 点 p 移 动 一 个 单位 到 吸收 壁 后 
动 。 


3.2 存在 吸收 壁 情况 下 的 有 限 限 步 一 维 随机 游 动 


问题 


» 


位 置 a 有 一 道 反射 壁 ， 求 点 P 经 过 nm 次 随机 游 动 出 现在 数 轴 


便 停 在 原 处 不 再 作 随机 游 


x 各 位 置 的 概率 。 


分 析 


经 过 n 次 随机 游 动 ， 点 P 出 现在 位 置 s 的 可 能 情况 〈 表 2 是 a=3 的 情况 ) 


S -10 -9 -8 = 人 -6 -6 -4 =8 = = 0 1 2 3 

1 1 1 

2 1 2 1 

3 1 3 3 1 

4 1 4 6 3 2 

5 1 5 10 9 7 

6 1 6 15 19 9 14 

7 1 7 21 34 28 37 

8 1 8 28 55 62 28 74 

9 1 9 36 83 117 90 176 

10 1 10 45 119 200 207 90 352 

表 4 
由 表 4 可 得 , 经 过 mn 次 随机 游 动 , 点 P 出 现在 位 置 s 的 可 能 情况 的 排列 满足 “ 缺 ” 的 “ 杨 
辉 三 角 ” 的 排列 。 

i . -10 9 -8 -=F -6 = -4 -3 =2 5 0 1 2 3 4 5 6 8 10 
1 

2 

3 

4 1 1 

5 1 5 中 

6 1 6 6 1 

WA 1 7 21 7 

8 L 8 28 28 8 1 

9 让 9 36 84 36 

10 1 10 45 120 120 45 10 1 

表 5 
对 比 表 4 和 表 5 可 得 , 表 4 是 由 表 1 减 去 表 5 而 得 , 并 且 对 应 吸收 壁 a 分 别 是 表 5 的 行 
和 。《【 见 表 5) 


不 妨 令 0<a<n，aEZ， 那 么 点 P 出 现在 位 置 s 的 概率 Ps 满足 


I 


s=2a-n 时 
P; (X=s) =P, (X=s) 


长 


当 2a-n 三 sa 时 
Ps (X=s) =Pi (X=s) Pi (X=2a-s) 
当 s=a 时 


Ps (X=s) 9 > Pi(X= s) SP (X=a) 


当 Ss 之 a 时 
Ps (X=s) =0 


我 们 可 以 把 吸收 壁 看 作 是 一 面 特殊 的 镜子 ， 那 么 点 P 出 现在 位 置 s 的 概率 就 等 于 点 P 出 现 
在 位 置 s 的 概率 与 位 置 s 关 于 吸收 壁 a 对 称 位 置 s” 的 概率 之 差 , 而 点 P 出 现在 位 置 a 的 概率 就 等 


于 点 P 出 现在 位 置 a 的 概率 再 加 上 两 倍 的 所 有 位 置 s 关 于 吸收 壁 a 对 称 位置 s” 的 概率 之 和 。 


点 P 出 现在 数 轴 x 各 位 置 概率 Ps 的 分 布 列 


当 n, a 同 为 奇偶 数 时 
—n —(n-2) … 2a-n-2 2a-n 2a-n+2 
(1 pn OE Ee ( n ] J 对 卫 ] . (三 je = ny n f 十 Pe 
1 a—1)? ? Ey S (9) (ja 
a-2 a 
nt+a 1 2 一 a 1 一 | | 站 n+a i n—a i 开征 沟 
jp jp? (jo? pp! 2. Pe pe es 
el | 0 | 
2 \ 2 \ 


当 n, a 不 同 为 奇偶 数 时 


—n =(n-2) … 2a-n-2 2a-n 2a-n+2 


(1 一 D) (rp ee | | Ea ( , ] 地 Go es ( 


a—1 


a-3 a-l 


n ‘ 加 2 n 区 a 2 n 2 
n+a-3|? (一 月 [ss (1—p) [sj (1 一 器 


D 十 a 一 3 2 一 a 十 3 mn 十 a 一 3 了 一 a 十 3 nt+a—3 n-at+3 nt+at+l et nd 


2 
‘1p 2 2 之 


全 四 
n+a++l Ee 
过 


2 


14 


3. 3 存在 反射 壁 情 况 下 的 无 限 步 一 维 随机 游 动 


当 x<a 时 


当 x=a 时 


当 x> 二 a 时 


3.4 小 结 


恨 据 1. 4 的 结论 和 3. 2 的 分 析 得 


fs (x)=f1 (x) -fi (2a-x) 


£0 72 A dx 


a 


fs (x) =0 


上 可 知 ， 所 得 存在 一 道 吸收 壁 的 情况 下 的 一 维 随机 游 动 


_ x2a)? 
“8 2 一 2 x<a 
_| V2x V2 
fs(x)= 
了 
?| ‘ec 2 dx -dirac(x 一 0) x>a 
V2T 
吻合 的 越 好 。 (dirac 即 狄 拉克 函数 ) 


图 6 


能 够 较 好 的 体现 


点 P 出 现 位 置 的 概率 密度 函数 


概率 分 布 ， 其 ! 


n 越 大 


结论 


经 过 多 种 情况 的 讨论 和 分 析 , 我 们 发 现 随 机 游 动 的 最 终 趋势 是 满足 正 态 分 布 的 。 而 在 反 
射 壁 的 干扰 下 正 态 分 布 会 在 反射 壁 附近 产生 聚集 现象 ( 见 图 7) ， 反 射 壁 离 初始 位 置 越 近 ， 聚 
集 现 象 越 明 显 。 在 吸收 壁 的 干扰 下 正 态 分 布 会 在 吸收 壁 附近 产生 排斥 现象 〈 见 图 8) ， 离 吸 
收 壁 越 近 以 及 吸收 壁 离 初始 位 置 越 近 ， 排 斥 现象 越 明 显 。 
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